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ABSTRACT 

For active galactic nuclei (AGNs) we study the role of the mechanism of quasi-linear diffusion 
(QLD) in producing the high energy emission in the MeV-GeV domains strongly connected with 
the submillimeter/infrared radiation. Considering the kinetic equation governing the stationary 
regime of the QLD we investigate the feedback of the diffusion on electrons. We show that this 
process leads to the distribution of particles by the pitch angles, implying that the synchrotron 
mechanism is no longer prevented by energy losses. Examining a reasonable interval of physical 
parameters, we show that it is possible to produce MeV-GeV 7-rays, strongly correlated with 
submillimeter/infrared bands. 

Subject headings: galaxies: active — gamma rays — submillimeter — infrared radiation 



1. Introduction 

According to the model of AGNs, cold material 
close to the central black hole forms an accretion 
disc, matter inside which due to the dissipative 
forces transports inwards causing the accretion 
disc to heat up. Such a hot material in turn can 
inverse -Compton scatter pho tons up to X-ray en- 



ergies (jBlandford et alJ ll990). From high energy 



astronomical sources a special interest deserve 
blazar type AGNs the standard model of which 
implies the presence of the supermassive black 
hole, surrounded by the accretion disc and ejecting 
twin relativistic jets. The observationally evident 
broadband emission spectrum of blazars is made of 
two components: the low energy (from radio to op- 
tical) domain attributed to the synchrotron emis- 
sion and the high energy (from A-rays to 7-rays) 
part formed by either the inv erse-Compton mecha- 
nism dBlandford et aljl990l) or the curvature radi- 
ation (|Gangadharalll996t iThomas fc Gangadharal 
20051 ). A recent investi gation of the p a rsec s cale 



jets is very important. Giroletti et al. ( 2010h 



gued that the high energy and radio emissions are 
strongly correlated. Our model on the other hand, 
as we will see, automatically provides a connection 



of radiation in high energy and radio domains. 
Magnetospheres of AGNs have strong magnetic 
fields, therefore, synchrotron cooling timescales 
are relatively short, leading to efficient energy 
losses. This in turn creates appropriate conditions 
for particles to transit to their ground Landau 
state. When this happens, relativistic electrons 
will move only along magnetic held lines without 
emitting in the synchrotron regime. 

Despite the very strong magnetic field, which 
prevents a continuous process of the synchrotron 
emission, there is a possibility to overcome the dis- 
sipativ e factors and maintain th e radiation mech- 
anism. iMachabeli fc Usov (Il979h have studied the 
cyclotron instability of two-component electron- 
positron plasma for the pulsar, NP 0532. It was 
found that the instability arises near the light 
cylinder surface (a hypothetical zone, where the 
linear velocity of rigid rotation equals exactly 
the speed of light) leading to a certain distri- 
bution of particles by pitch an gles and the con- 



sequen t synchrotron radiation. lLominadze et al 
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(Il979h considered the magnetospheres of the pul- 
sar NP 0532 and the Crab nebula, studying the 
generation of waves from optical to gamma-ray 
dom ains. The similar appro ach was presented 
by iMalov fc Machabelil (|200lh where the QLD 
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was applied to the radio pulsars. The authors 
found that the transverse momenta of relativis- 
tic particles induced by the cyclotron instability 
caused the stable non-zero pitch angle distribu- 
tion maintained by means of the QLD. Analyz- 
ing the data obtained from MAGIC Chercnkov 
teles cope between 20 07 October and 2008 Febru- 
ary (jAliu et al.ll2008l ) we found that, the observed 
coincidence of signals in the optical and 7-ray 
domains are easily explained by the QLD pro- 
cess, which leads to the increase of the pitch 
angles, making the synchrotron p rocess feasible 
(|Machabeli k Osmanovll2009l l2010( ). 

In the magnetospheres of AGNs the magnetic 
fields are of the order of 10 4 G (Thorne et al 



condition 



19881 ). close to the supermassive black hole, and 
100G — 300G close to the light cylinder sur- 
face. Therefore, the aforementioned QLD mech- 
anism could be of great importance for AGNs 
as well. For this purpose by considering the cy- 
clo tron instability excited in th e radio domain, 
in ( Osmanov k Machabeli 2010l) we studied the 
quasi-linear interaction of proper modes of AGN 
magnetospheric plasmas with the resonant plasma 
particles investigating the QLD in the context of 
producing the soft and hard X-ray emission from 
AGNs. Under favorable conditions this mech- 
anism could also be efficient for explaining the 
MeV-GeV energy synchrotron emission, strongly 
connected either with the submilimeter radio 
band, or with the infrared emission induced by 
the cyclotron instability. This will be the subject 
of the present paper, which is organized as fol- 
lows. In Section 2 we describe our model, in Sect. 
3 we apply the mechanism of QLD to AGNs and 
in Sect. 4 we summarize our results. 

2. Main consideration 

In general, AGN magnetospheres consist of rel- 
atively low energy particles and very high energy 
particles (electrons). Therefore, in the framework 
of the model we consider the plasma composed of 
two components: a) the so-called plasma compo- 
nent with the Lorentz factor, j p and b) the beam 
component with the Lorentz factor . 7^ (7?, ^> 7 n 



Such a system as was shown by (jKazbegi et al 



1992) undergoes the cyclotron instability induced 
by the Doppler effect with the following resonance 
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0, 



(1) 



where k., is the longitudinal (parallel to the back- 
ground magnetic field) component of the wave vec- 
tor, k x is the component along the drift, V.. is the 
longitudinal component of plasma flow velocity, 
u x = cV 'Yb/ P&B is the drift velocity of particles, 
c is the speed of light, p is field line's curvature 
radius, tug = eB/mc is the cyclotron frequency, B 
is the magnetic induction and e and m are elec- 
tron's charge and the rest mass, respectively. Pos- 
itive sign corresponds to the damping of the ex- 
cited modes, whereas the negative sign relates to 
the unstable mode. One can show that, when the 
aforementioned resonance takes place, the trans- 
verse waves with the dispersion relation 



fee (1-5), 5 = 



(2) 



are induced, k is the modulus of the wave vector, 



y / 47rn p e 2 /m is the plasma frequency and 



is the plasma density. From Eqs. (|1I2|) is clear 
that the excited cyclotron frequency is given by 
(jMalov k Machabeli|[200ll) 



U! 



MB 
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(3) 



In spite of the resonant character of the cyclotron 
modes, the corresponding frequency is not well 
peeked, because ui depends also on the Lorentz 
factors of the resonant (beam) particles, that do 
not have narrow energy spectra. It is worth not- 
ing that unlike the synchrotron mechanism (A < 

— 1/3 

n p ), (A is the wavelength), where radiation pro- 
cess can be described by a single particle approach, 
excitation of the aforementioned waves is a collec- 

— 1/3 

tive phenomenon (A > n p ) , which in its turn is 
a direct consequence of one-dimensionality of the 
distribution function. On the other hand, such 
a behaviour of this function is guaranteed by the 
strong magnetic field that forces particles to move 
along the field lines. 

By the first sentence after Eq. 3 we would like 
to note that despite the resonant character of the 
process, the corresponding frequency is not well 
peeked, because the frequency depends also on the 
Lorentz factor of the resonant particles, and the 
corresponding interval of 7&-S is not narrow 
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In general, two dissipative factors lead to a de- 
crease of the pitch angle. The force that provides 
conservation of adiabatic invariant I = 3cp 2 /2eB 



particles and 



in no n-uniform magnetic field ([Landau fe Lifshitz 
197lh [see also (|Osmanov fc Machabehll2010l) ] 



ttic 

G_L = IbW, 

P 



G x = — 7^ 2 (4) 
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and the radiation reaction force ( Landau fc Lifshitz 



19711) 



F x = -aV(l + 7 h V), ^ = 



-ocilV, (5) 



where a = 



2e 2 c4/(3c 2 ) 



and "0 is the pitch an- 



gle. Only under the action of these forces the 
pitch angles tend to decrease, inevitably killing 
the synchrotron emission. In reality the situation 
is principally different, because the excited rela- 
tively low frequency waves (in our case submillime- 
ter/infrared) by means of the cyclotron resonance 
leads to the QLD. Unlike the dissipative effects 
of (F, G), diffusion creates non-zero distribution 
function with respect to pitch angles, and supports 
the synchrotron emission. It is clear that under 
favorable conditions the effect of quasi-linear dif- 
fusion may balance the dissipation, therefore, our 
objective is to find the distribution of particles by 
pitch angles, estimate their mean value and an- 
alyze the corresponding synchrotron emission en- 
ergy. On the other hand, since the QLD results 
from the feedback of the cyclotron modes, apart 
from the high energy emission, the system will also 
be characterized by the low energy radiation. 

In ( Osmanov fc Machabelil 201dh we studied 
the role of the QLD in producing the X-ray emis- 
sion by means of ultra-relativistic electrons in 
AGN magnetospheric flows. It was shown that 
the cyclotron resonance provides emission in a low 
energy domain - radio band. Unlike the physi- 
cal conditions (\G i I > |F| I and \G„ | <C \F {1 \) 
considered in ( Osmanov fe Machabelil l201fj) , in 



the present paper we examine physically differ- 
ent regime, |Gj_| < \F±\ and |GJ < |FJ, which 
reduces the stationary kineti c equation governin; 
the q uasi-linear diffusion to ( Malov fc Machabel 
20011) 
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dtp 
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where / = f(ip) is the distribution function of 
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is the diffusion coefficient. \Ek\ 2 is the energy 
density per unit wavelength, therefore the energy 
density of the cyclotron waves are of the order 
of \Ek\ 2 k. We assume that ~ 50% of the reso- 
nant plasma energy, roc 2 n&7& , converts to waves 
(jOsmanov fc Machabelil I201C), then for an ex- 
press ion of \Ek\ 2 one obtains (jMalov fc Machabeli 
20011) 



\Ei.\ 2 



mc nb"fb 
2w ' 



(8) 



The distribution function obtained from Eq. ([7|) 
is given by 

f(fl>)=Ce- A *\ (9) 



where 



.4 



4D, 



C = const. 



(10) 



UnUk e the work presented in ( Osmanov fc Machabelil 
2010), due to the different regime, f(ip) behaves 
as e~ A ^ instead of e~ Al ^ [see Eq. (11) in 
(jOsmanov fc Machabeh1l2010l )]. 

As we see, particles are distributed by the pitch 
angles, therefore, the electrons will emit via the 
synchrotron process without damping. 

3. Discussion 



In this section we apply the model of the QLD 
to the light cylinder lengthscales of typical AGNs. 

For explaining the high energy radiation, it is 
strongly believed that AGN magnetospheres are 
consist of highly relativistic electrons. This fact 
sets another problem - how to accelerate these 
particles to such high energies? In general, there 
are several mechanisms, which may account for 
the efficient acceleration of electrons. Indeed, as 
is shown in a seri es of papers, Fermi-typ e ac- 
celeration process ( Catanese fc Weeks 1999;), re- 
acceleration of el ectron-positron pairs as a feed- 
back mechanism ( Ghisellini et al. 19931) and cen- 
trifugal acceleration dMachabeli fc Rogava] I l994 ; 
Osmanov et al. 120071 : iRieger fc Aharonianl 12008 ; 
OsmanovT 2010 ) may provide very high Lorentz 
factors of the order of jb ~ 10 5-9 . Therefore, 
in the framework of the paper an existence of 
such particles is assumed to be as a given fact. If 



3 




Fig. 1. — Behaviour of CMeV versus jb for dif- 
ferent values of n&. The set of parameters is 
L = 10 45 er 5 /s, tt = 3 x lO^rad/s, 7p = 200 
and rib = {5; 10; 15}cm~ 3 . As is seen from the 
plots, relativistic electrons with Lorentz factors 
76 = {1 — 2} x 10 s may provide the high energy 
radiation in the MeV-GeV domain. 
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Fig. 2. — Behaviour of €MeV versus w for dif- 
ferent values of rib- The set of parameters is 
L = 10 45 er 5 /s, S] = 3x lO^rad/s, 7p = 200 
and rib — {5; 10; 15}cm -3 . We see that the high 
energy radiation is strongly connected with the 
submillimctcr/low infrared emission. 



we suppose an isotropic distribution of relativistic 
electrons, one can estimate the synchrotron cool- 



ing timescale (jOsmanov fc Machabelill2010l ) 
£ cooi « 5 x HP 3 x f^Vx (— )s. (11) 



B 



The value of the magnetic in duction is given by 
(jOsmanov fc Machabeh1l2010[ ) 



260 x 



L 



(l0 45 er ff /i 



1/2 



3 x 10-V 



G, 
(12) 

where L is the bolometric luminosity of the AGN 
and ri c — c/£l is the light cylinder radius (a 
hypothetical zone, where the linear velocity of 
rigid rotation exactly equals the speed of light). 
f2 is the magnetic field lines' angular velocity of 
rotation, normalized to the value 3 x 10~ 5 s~ 



( Belvedere et al. 1989). Throughout the paper we 



assume that the magnetic field is robust enough 
to maintain the frozen-in condition in the magne- 
tosphere of the AGN. Indeed, as one can see, for 
the typical magnetospheric parameters, 7f, ~ 10 s , 
rib ~ 10crn~ 3 , the following condition Bf c /8n > 
7f,mnf,c 2 is satisfied, which means that the plasma 
particles will be forced to follow the rigidly ro- 
tating field lines. During such a motion, espe- 
cially on the light cylinder lengthscales, the elec- 
trons will undergo the centrifugal force accelerat- 
i ng them up to very high Lorentz factors ~ 10 8 ~ 9 
(jOsmanov et al.ll2007t iRieger fc Aharonianll2008j) 
It is clear from Eq. (fTTj) that for a certain 
class of physical parameters the synchrotron cool- 
ing timescale is of the order of 5 x 10 _4 s. On the 
other hand, the kinematic timescale of the sys- 
tem, tkin ~ ric/c ~ 3 x 10 4 s is by many orders of 
magnitude bigger than t coo i, which in turn means 
that without the quasi-linear diffusion, particles 
very soon would stop emitting in the synchrotron 
regime, after transiting to their ground Landau 
level. 



Kazbegi et al 



(|1992f) showed that the anoma- 
lous Doppler effect generates the cyclotron wave s 
with the frequency ( Osmanov fc Machabeli 2010f) 



w w 6.8 x 10 9 x 
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rib 



(13) 
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leading to the process of the quasi-linear diffusion, 
which, despite the efficient dissipative factors, cre- 
ates the pitch angles. 

To demonstrate the present model, we con- 
sider an AGN with the bolometric luminosity 
L = lO^erg/s. Let us examine the following 
parameters SI — 3 x 10~ 5 rad/ s, jb — fO 8 and 
rtf, = 10cto~ 3 . Since the particles are distributed 
by the pitch angles [see Eq. ©], for analyzing the 
synchrotron emission it is reasonable to estimate 
a mean value of ip> 



aforementioned condition reduces to 



0.5 



Trr- ( 14 ) 



Then one can show from Eq. (114)) that for the 
aforementioned parameters the pitch angle is of 
the order of 8 x 10~ 3 rad, therefore, relativistic 
elect rons will inevitably emit p hotons with ener- 
gies ( Rvbicki fc Lightman 19791) 



(15) 



e eV w 1.2 X 10" 8 B7 2 sinV'. 



After substituting the value of ip in Eq. (|15p . we 
see that the synchrotron emission generates radi- 
ation in the MeV domain. 

The quasi-linear diffusion works if the cyclotron 
modes are excited, therefore, it is essential to es- 
timate the timescale of the corresponding insta- 
bility (ti ns ) and compare it with the kinematic 
ti mescale of the system . According to the work 
of iKazbegi et al.l (1992) the growth rate of the in- 
stability is given by 



and 



r = 7T 



2oj7 p 5 ■ 



if 



if 



luj 
2 c 2 



«<5 



1 u l r 

2 cr 



(16) 



(17) 



where Wb = \J Aim^e 2 / 'm is the plasma frequency 
of beam electrons. It is easy to show that for 
n b = 10cm- 3 , 7 P = 200 (see Fig. HJ, Vj, - c 
and p ~ R g , one obtains u 2 /(2c 2 ) -C 1, implying 
that the increment of the instability is given by 
Eq. (1161) . The cyclotron resonance makes sense 
if Uns/tkin < 1) then, by taking into account the 
definition of the kinematic timescale, ri c /c and the 
instability timescale, 1/r, one can show that the 



3.5x10" 



xfl^l 1 

V rib 



V100/ \7b 

(18) 

As is clear from Eq. (TT5|) . the condition is very 
sensitive to the Lorentz factor of the plasma com- 
ponents, and for relatively higher values of 7 P the 
condition will violate. The upper limit of j p , when 
the condition is still valid is of the order of 300 for 
7i, ~ 10 8 and rib ~ IOcto -3 . 

For studying the efficiency of the QLD, we ex- 
amine the following set of the parameters L = 
10 i5 erg/s, Q = 3 x 10" 5 rad/s, 7p = 200 and 
rib = {5; 10; 15}cto~ 3 and the results are demon- 
strated in FigQ] where we show the behaviour of 
e Afev versus 7b. From the plots is clear that e MeV 
is a continuously increasing function of the beam 
Lorentz factor, which is a natural result of the fact 
that more energetic particles produce more ener- 
getic photons. The behaviour of e McV versus rib 
is different, more dense beam electrons produce 
photons with lower energies. This can be seen 
from Eqs. (|10I14|) : tp ~ y/D ±± , which by com- 



bining with D ±± 



[see Eqs. (|7I8[) ] confirms 



the dependence e MeV (n.b). According to the re- 
sults demonstrated in the figure, relativistic elec- 
trons with Lorentz factors 7b = {1 — 2} x 10 s may 
provide the high energy radiation in the MeV-GcV 
domain. 

Since the generation of the synchrotron emis- 
sion strongly depends on the cyclotron waves, we 
also investigate the behaviour of e MeV versus u). 
Figure [2] shows the function e MeV (ui) for several 
values of rib. The set of parameters is the same 
as in the previous figure. As is clear from the 
plots, the {200 - 1200}MeV radiation is strongly 
connected with the submillimeter (~ [0.3 — 3] x 
10 12 iJz) and low infrared (~ [3 - 3.8] x 10 12 Hz) 
emission. 

Another important parameter, the physical sys- 
tem depends on, is the Lorentz factor of the 
plasma component. Therefore, it is reasonable to 
demonstrate the function e M(sV (w) for different val- 
ues of 7 p . These results are shown in Fig. [3l 
where the set of parameters is L — 10 erg/ s, 
n = 3 x lQ- 5 rad/s, 7 P = {200; 250; 300} and 
rib = {5; 10; 15}cm" 3 . As we see, the excited in- 
frared domain extends up to ~ 10 13 Hz, which 
is strongly connected with high energy emission 



5 



(lOOMeV). As is clear from the figure, higher val- 
ues of 7 P correspond to lower synchrotron ener- 
gies. Indeed, by taking into account the relation 
rj) ~ y/D ±± combined with Eqs. (|7I8| one can see 
that ijj ~ I/72. 

As is clear from the results, the quasi-linear 
diffusion together with the cyclotron instability 
may guarantee production of high energy radia- 
tion in the MeV-GeV domains strongly connected 
with submillimeter/infrared emission. The ma- 
jor diffe rence in results from our prev ious work is 
that in (jOsmanov k Machabelill2010l ) we studied 
physical conditions leading to excitation of X- 
rays connected with the relatively low frequency 
radio band (KHz-MHz), whereas in the present pa- 
per both energies (produced as by synchrotron as 
by cyclotron mechanisms respectively) are much 
higher. This investigation sets another problem: 
since the QLD is a feasible mechanism provid- 
ing the aforementioned high energies, one of the 
important next steps could be testing of MeV- 
GeV AGNs exhibiting an efficient submillime- 
ter/infrared radiation and see if the strong cor- 
relation is observationally evident. This in turn, 
could be a certain test for estimating the AGN 
magnetospheric parameters, such as the density 
and the Lorentz factors of plasma component and 
beam electrons. A particular future objective is to 
investigate theoretically the radiative signatures of 
both high and low energy emissions respectively. 
This work we are going to perform sooner or later. 

4. Summary 

The main aspects of the present work can be 
summarized as follows: 

1. Mechanisms producing strongly connected 
high and low energy radiation was stud- 
ied by taking into account the QLD in 
the AGN magnetospheres. Considering 
a physical regime different from that of 
( Osmanov k Machabelil l2010h . we investi- 
gate the efficiency of the QLD in a region 
close to the light cylinder surface. 

2. For the considered physical parameters, it 
has been shown that the cyclotron insta- 
bility appears for relatively low frequency 
range, producing radiation in the submil- 
limeter/infrared domains. On the other 



hand, despite the short cooling timescales, 
the effect of diffusion on particles recreates 
the pitch angles and produces the high en- 
ergy radiation in the MeV-GeV bands. 

3. The problem was studied versus three ma- 
jor magnetospheric parameters: the beam 
and the plasma component's Lorentz fac- 
tors and the beam electrons' density. It 
was shown that the photon energy, e MeV , 
is a continuously increasing function of the 
beam Lorentz factor and the beam density. 
Contrary to this, by increasing the plasma 
component Lorentz factor, the correspond- 
ing photon energy decreases. 
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Fig. 3. — Behaviour of eueV versus co for dif- 
ferent values of j p . The set of parameters is 
L = 10 45 er 3 /s, il = 3 x lO^rad/s, lp = 
{200; 250; 300} and n b = {5; 10; 15}cm" 3 . The 
results show the strong connection of {100 — 
800}MeV radiation with the emission in the fre- 
quency range {0.05 - 1} x 10 13 Hz. 
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